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Abstract: We propose a scheme for efficient cavity-enhanced nonlinear 
THz generation via difference-frequency generation (DFG) processes 
using a triply resonant system based on photonic crystal cavities. We show 
that high nonlinear overlap can be achieved by coupling a THz cavity 
to a doubly-resonant, dual-polarization near-infrared (e.g. telecom band) 
photonic-crystal nanobeam cavity, allowing the mixing of three mutually 
orthogonal fundamental cavity modes through a x^^'' nonlinearity. We 
demonstrate through coupled-mode theory that complete depletion of the 
pump frequency - i.e., quantum-limited conversion - is possible in an 
experimentally feasible geometry, with the operating output power at the 
point of optimal total conversion efficiency adjustable by varying the mode 
quaUty (g) factors. 
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1. Introduction 

Nonlinear optical frequency conversion is widely used for the generation of light in parts of the 
spectrum for which there are no convenient sources iHl . In particular, nonlinear processes are re- 
garded as a promising route to generation of coherent radiation in the terahertz (THz) frequency 
range |l2]|3]|4]|5]|6]|2l[8]|9][T0l. Recently, there has been renewed interest in cavity-enhanced 
nonlinear frequency conversion, as improved designs and fabrication techniques have paved 
the way for the realization of wavelength-scale cavities, thus allowing efficient conversion at 
increasingly low powers ||2[3]|l|5]|6l[IIl[I2[I3IIlll5l[l6l[l7l[l8l. 

In this paper, we propose a scheme for efficient THz generation based on second-order j^^' 
difference frequency process in a triply resonant structure consisting of two nested photonic- 
crystal cavities, as shown schematically in Fig. 1. The first cavity, designed to operate at near- 
infrared (NIR) wavelengths and realized in a nonlinear material (e.g. III-V semiconductors), 
produces nonlinear-polarization at THz frequencies by difference-frequency generation (DFG). 
The cavity is placed in the proximity of a second, much larger, single-mode cavity with the res- 
onance in the THz range. We utilize the fact that the THz wavelength. At (~ 150/im), is on 
a vastly different length scale from the pump and idler wavelengths, Ai and A2 respectively 
1.5/im), allowing us to effectively decouple the cavity designs. Our system can be viewed 
as a classical version of a two-level atom strongly coupled to the optical cavity. In our design. 




Fig. 1. Design of tlie triply resonant system of coupled photonic crystal nanobeam cavities 
for efficient THz generation. A dual-mode PhC nanobeam cavity is suspended just above 
tlie THz cavity (there is a 1/im air gap in our simulations) near its field maximum, (a) 
Normalized mode profile (E. j) of the THz mode. Dimensions are shown as a function 
of the THz wavelength. Dimensions are shown here for GaAs, but in principle any mate- 
rial and design can be used as long as it has a THz mode with the correct polarization, 
(b) Normalized nonlinear mode-product (£v,TE^.r.TM)of the TE-like and TM-like modes 
in the dual-mode cavity which overlaps with the THz mode shown. This product has a 
single sign, indicating constructive overlap of the two modes. Dimensions other than the 
beam thickness are shown as a function of the lattice period a. (c): Nanobeam thickness 
(left) and Bragg mirror periodicity, (right) for the NIR PhCNC plotted as a function of the 
THz difference frequency. The TE-like mode frequency is fixed at 200THz (Aje = 1.5/im, 
Ate < ■^tm)- 

we suspend the telecom band cavity above the THz cavity near the field maximum |fT9l . suf- 
ficiently far so that the telecom modes are not affected by the presence of the THz cavity, but 
sufficiently close that the long-wavelength THz mode extends substantially into this cavity (e.g. 
1/im above). In this way we can achieve a high overlap between the three fundamental cavity 
modes by mixing three mutually orthogonal polarizations through the strong Xij/^ (' 7^ j ¥^ k) 
terms in the nonlinear susceptibility of III-V semiconductors (e.g. GaAs, GaP) ||2]|2ll- This 
scheme takes advantage of our recently demonstrated doubly-resonant ultrahigh-Q photonic- 
crystal nanobeam cavities (PhCNC) with mutually orthogonally polarized modes (i.e. one is 
TE-like one is TM-like) ll9l 1201 . Both cavities are preferentially coupled to a corresponding 
waveguide extending from one end, with the length of the Bragg mirror of holes at one end 
used to tune the strength of the coupling. The mirror at the other end is made sufficiently long 
such that leakage in that direction is much smaller than out-of-plane losses. Since the relevant 
nonlinear material is in the NIR cavity, the THz cavity can be composed of any material, a 
useful degree of freedom considering the scarcity of low-loss THz materials. 



2. Model 



We employ the coupled-mode theory (CMT) framework described in Ref. lfT4l to analyze triply 
resonant THz generation by DFG in a wavelength-scale cavity with a non-resonant j'^' non- 
linear susceptibility. Our cavity has three resonant modes at frequencies ©i (pump), ©2 (idler), 
and a>r (THz signal), which satisfy ©r = ©1 — ©2 and each mode is coupled to a waveguide. 
As described in Ref. 1141 . the conversion efficiency of the system is determined by the input 
powers, {Pi,P2), the cavity mode quality factors, Qk{k~ 1,2,T), the cavity waveguide coupling 
strengths F^. (defined as F^. = yi^ J where 7/. j is the leakage rate into the waveguide and ji^ is 
the total leakage rate including linear losses), and the nonlinear coupling constant j3, defined in 
Ref. lfT4l . As discussed in Ref. lfT3l . avoiding mode symmetries and polarizations that lead to 
zero coupling constant, as well as optimizing the value of j3, is analogous to phase-matching in 
the nonlinear mixing of propagating modes. For our type of design, since the THz field strength 
is nearly constant across the much smaller NIR cavity, j8 is well approximated by: 



TM| 



where Eje, Etm- ™d Ex are the field profiles of the cavity modes, denotes spatial integration 
only over the nonlinear dielectric, and is a dimensionless constant that quantifies the THz 

cavity contribution to the nonlinear overlap, defined as: K-y = A;^^^£'-.t,nir/[\/ / tj^^rerlExp], 
where S^^t.nir is the constant THz field amplitude inside of the NIR cavity. 

As described in Ref. |fT4l . quantum conversion efficiency is quantified by the parameter: 

eQ = -^^. (2) 
*^ ©tFiFt Pi 

This parameter describes the number of THz photons collected per photon of signal (©1) input 
that couples into the cavity (factoring in losses). As expected, this parameter is bounded in 
the steady state by < 1 lfT4l . Since the power efficiency (Poutj/Pi) is scaled by a factor 
of 1 /FjFx, a decrease in FiF-p causes a decrease in the maximum attainable power efficiency 
(always occurring at E^ = 1). Therefore over-coupling of the cavity modes to the corresponding 
input/output waveguides is imperative for optimal performance lfT4l . As derived in Ref. lfT4l . 
the efficiency of conversion for any triply resonant DFG system can be expressed as a function 
only of the input powers, normalized against the critical powers (Pi.crft, ^'2,crit) defined as, 

where Q = QiQiQt- This means that only the product of the three Q factors contributes 
to efficiency and the distribution of this product is unimportant. Figure 2a shows the stable 
quantum efficiencies reached after a simple continuous-wave (CW) step excitation as a func- 
tion of normalized input powers. The dotted black line denotes the onset of bistability to its 
right lfT4l . Quantum-limited conversion {E^ = 1 ) can be achieved when the input powers satisfy 
Pi/ P2,cni = (1 — A/4f'i,crit)^- We Can also define a total-efficiency parameter that is proportional 
to the ratio of THz output power to total input power at both telecom frequencies: 

* ©rFirTA+^'2 Pi+Pi' 
Quantum efficiency and total efficiency are equivalent in the asymptotic limit, f 2 ^ (P2 ^ 
^2,crit> P2 ^ P\)- In this limit, the conversion efficiency is mono-stable lfT4l . Total power con- 
version is optimized when P2 <C ^'2,crit- P\ = 4/"! .crit- One can calculate a trajectory in the (Pi ,P2) 



space that maximizes Eg. This corresponds to the optimal operating conditions for a single 
cavity (fixed Q, {(0^}, j3) operating through a range of powers. The white dashed curve in 
Fig. 2a shows the ideal operating conditions for CW THz generation (the limit coi w ©2). For 
Pi > 4Pi.crit, there do exist solutions with slightly higher efficiency than on this trajectory, but 
these solutions are only stable for certain g-f^ctor combinations and cannot be excited with a 
step CW excitation. Keeping P2 <<Pi for all Pi > 4Pi crit ensures mono-stable behavior in this 
region lfT4l . 

3. Design 

In order to estimate the coupling constant, j3, we consider a geometry similar to that shown in 
Fig. 1: a PhCNC, with closely-spaced TE-like and TM-like modes in the near infrared (NIR) 
telecom spectral region {(0\/2n, (ol/ln ^ 200THz) ||20:| positioned 1/im above a larger THz 
cavity. Our cavities, both telecom and THz, are based on a free-standing ridge waveguide pat- 
terned with a one-dimensional (ID) lattice of holes, utilizing a tapered photonic-crystal cavity 
approach that we and others have recently demonstrated ||9l |22l |23] |24l |25l . We consider the 
case where both cavities are composed of GaAs (hnir ~ 3.45, hthz ~ 3.6) ifTOl 1211 . and the 
telecom cavity is suspended slightly above the THz cavity at a height, d = 1/xm <C Xj. By en- 
gineering the "impedance matching" between the cavity mode and evanescent Bloch mode that 
exists in the photonic crystal (Bragg) mirror we have demonstrated, both experimentally and 
theoretically, Q factors on the order of 10^ |l9l|23l. A 3D finite-difference time-domain code 
(3D-FDTD) 1261 was used to design independent telecom and THz cavities 123, and a CMT 
model (as described above) was used to describe the nonlinear frequency conversion. 

We use the design algorithm detailed in Ref. |9!] for the design of an NIR PhCNC support- 
ing orthogonal TE-like and TM-like modes, which is necessary to achieve efficient nonlinear 
conversion as dictated by the x'"^^ tensor of a III-V semiconductor material. We fix the ridge 
width equal to the regular lattice period (a) and use circular holes in the Bragg-mirror region 
with radius, R = 0.3a. The period and hole radius are then tapered linearly in reciprocal space 
to 0.84 times the normal value over 8 periods ||9|- All cavity dimensions are fixed with respect 
to the Bragg-mirror lattice period (a) except for the ridge thickness (or height), which is used to 
tune the frequency spacing between the TE-like and TM-like modes (see Fig. Id and Ref. ||9l). 
Here we vary the cavity dimensions so that the TE-like resonance is always at 200THz (1.5/im) 
and the TM-like resonance is at a slightly lower frequency, giving the desired THz difference 
frequency. The TM-like resonance can also be tuned above the TE-like resonance in frequency, 
but we find that this gives poorer TM Q factors and lower values of j3 as the TM-like bandgap 
shrinks in this direction ||9l- The NIR cavity parameters are plotted in Fig. Ic as a function of 
the desired THz frequency. 

For the THz mode, we consider a TE-like mode of a second GaAs PhCNC, based on the 
design shown in Fig. la, having elliptical holes (7?a = 1.57?b) and a four-hole, linear taper of 
the lattice period and hole diameters as was done with the NIR cavity. Fig. la shows the cavity 
dimensions as a function of the free-space wavelength of the resonance. This cavity can be 
scaled to fit the desired THz frequency and has Kt = 0.65 for any scaling. This cavity could in 
principle be replaced by any type of THz cavity having the same value of Kj, and all operating 
properties described below would be unchanged. This degree of freedom is a key advantage of 
our decoupled-design paradigm. The power at which maximum conversion efficiency occurs 
(normalized Pi = 4Pi crft) for a given THz frequency depends on the product of the three Q 
factors and the overlap, j3. Fig. 2b shows the dependence on the THz frequency in our design 
scheme of j3 and the Q-product required to produce ImW of THz radiation with optimal effi- 
ciency. Figure 2c shows the optimized NIR input power and optimized THz output power as a 
function of the g-product for our system for fj = 2THz, j3 = 3.5J^'/^. 




Fig. 2. (a) Step-excitable quantum efficiency (E^) for stable CW THz generation, plot- 
ted as a function of powers of the pump and idler normalized against the critical powers 
Pk.cnt (see Eq. l[3)). The solid line denotes the critical relationship between input powers 
where = 1 is possible (Pi = 4/"! cit). The black dotted line denotes the onset of bista- 
bility. An inset showing the other stable solution is shown in the top right corner 1141 . The 
white dashed line shows the optimal operating conditions for maximum total conversion 
efficiency (b): Performance parameters of our nested PhCNC design (GaAs) as a function 
of the THz resonance frequency: nonlinear overlap, j8 (left), and the Q-factor product re- 
quired for ImW of THz power to be generated from a pump at 200THz with quantum 
limited efficiency (right), (c): Dependence of the input power (solid line) yielding optimal 
efficiency (Pi = '^Pi^rit) and the corresponding THz output power (dashed line) on the cav- 
ity g-factor product (Q), for coupled GaAs THz and dual mode NIR PhCNCs (C0i/2k, 
0>2/2k ~ 200THZ, 0>y/27Z ~ 2.0THz, jS ~ 3.57"'/^). (d): THz output power as a function 
of input power in this geometry for, Q ~ 2.5xl0'^. The dotted line shows the quantum 
limit. 



The Q factors in our cavities were originally designed to be limited by out-of-plane scatte- 
ring II9I I22II . Overcoupling to the input/output waveguide is then achieved by reducing the length 
of the Bragg mirror on one side. This reduces the maximum values of Q by an order of mag- 
nitude, but minimizes the reduction in efficiency caused by losses 1 14], and enables operation 
close to quantum-limited conversion. We find that this variation of the Bragg mirror lengths has 
negligible effect on the overlap j3 for an order of magnitude variation in the Q-factors. For some 
applications where a single cavity must operate across a range of different power levels, the best 
performance is achieved by varying the two input powers to fit the trajectory described in the 
previous section (see Fig. 2a). For example, this optimal output power is shown in Fig. 2d for 
a cavity optimized to produce ImW of power at 2 THz. This optimization requires a g-factor 
product of Q ~ 2.5 x lO'^. When limited by out-of-plane losses, the Q factors of our TE/TM 
dual-mode NIR cavity are Qi (TE) = 2 x 10^, ^(TM) = 3 x 10*^ for = 2THz. Reduction of 
each mode's Q factor by an order of magnitude to overcouple the waveguide requires a THz Q 
of only ~ 300. This low requirement is important since realistic THz cavities will be limited 
by material losses in the THz range. While GaAs has a low absorption coefficient compared 
to most materials in the THz range, its linear loss rate still corresponds to a g factor of about 
1 .5 X 10^ lITOl . which is well below the scattering-Umited Q of our design (Qx.scat = 1 -4 x 10^). 
THz losses are likely to put an upper bound on any THz cavity Q factor in practice, and further 
reduction of this limit is then required for overcoupling. However, a large advantage of this 
geometry is that all nonlinear processes take place outside of the THz cavity material in the 
NIR cavity, so that any design and any material can be used for the THz cavity, with the per- 
formance of any THz cavity coupled to our dual-mode nanobeam cavity purely determined by 
the cavity's Kj value. Thus, this loss limit can be increased with the discovery of more efficient 
materials in the THz wavelengths. 

4. Arrays 

Considering that high-Q NIR modes are needed to completely convert even relatively high 
pump powers (see Fig. 2), the material damage threshold is likely to become an important 
design parameter The vast difference in length-scales in our design can be further exploited 
to circumvent this problem. This is because the NIR cavity is so small compared to the THz 
mode profile that many NIR cavities can be fabricated on top of a single THz cavity, all within 
the fundamental anti-node. In this limit, where the NIR cavity occupies a very small volume 
of the THz mode and where there is good overlap between the NIR modes, the NIR mode 
volumes have a negligible effect on j3 [see Eq. ^]. Therefore, we can consider an array of 

NIR cavities on the anti-node of the THz cavity (whose inputs are phase-locked), shown in 
Fig. 3a, as a single NIR cavity. This cavity has the same Q factor and a times higher mode 
volume, but will have the same j3 /kj, yielding an unchanged value of j3 provided that the array 
still fits entirely inside a region of the fundamental anti-node with roughly constant THz field 
amplitude. Thus, this cavity will have the same optimal THz output power, but with a damage 
threshold that is increased A^-fold. 

Quantum-limited efficient conversion can be achieved over a broad range of powers from 
a single chip by integrating many triply resonant devices as shown in Fig. 3b. The total size 
of a single THz device is given by the size of the THz cavity, which can be sub-millimeter 
(i.e., on the THz wavelength scale). Thus, hundreds of these devices can be integrated on a 
single square-centimeter chip. Large arrays allow not only for high-power efficient generation, 
but also maximally efficient generation over a broadly tunable range compared to what can be 
achieved from a single device. To see this, suppose a single device emits Poutj of THz power 
with E"^"^ = 1 . An array of of these devices therefore can emit A^Pout.x with = 1 , but can 
also emit at this efficiency over the whole range Pomj > P ^ NPomj by varying the number 




Fig. 3. (a): Large phased array of NIR cavities placed on a single THz cavity. This 
array will have the same THz- generation performance properties as the corresponding 
device having only a single NIR cavity, however the damage threshold input power will 
increase by a factor of N. (b): High-power tunable generation is possible by incorporating 
many (> 100) THz devices on a single chip. An array of N devices that emit Pout.T most 
efficiently allows quantum limited conversion over the full range of output powers from 
^out.T to A'-Pout.T- The specific scale shown in both (a) and (b) is for Aj = 150/im. 

of devices in the array that are used. 

5. Conclusion 

We have proposed an experimentally feasible platform for high efficiency THz generation in 
triply resonant photonic crystal cavities. By placing a dual-mode photonic-crystal nanobeam 
cavity with closely spaced TE-like and TM-like resonant modes in the telecom band near the 
field maximum of a much larger THz cavity mode, we have shown that a high nonlinear overlap 
of mutually orthogonal modes can be attained in many nonlinear materials (eg III-V semicon- 
ductors). Total power conversion approaches the quantum limit at a critical power which can 
be tuned by varying the product of the three mode Q factors. Coupling of many NIR cavities 
to a single THz cavity increases the material damage threshold without changing the device 
performance properties. High-efficiency power output can be further increased and tuned by 
fabricating large arrays of these triply resonant devices on a single chip. 

6. Acknowledgements 

We would like to thank J.D. Joannopoulos and Marin Soljacic for helpful discussions. IBB and 
MWM wish to acknowledge NSERC (Canada) for support from PGS and PDF programs. This 
work is supported through NSEC at Harvard, by the Army Research Office through the ISN 
under Contract No. W91 1NF-07-D-0004, and by US DOE Grant No. DE-FG02-97ER25308 
(ARW). 



